JET PRODUCTION IN ep COLLISIONS 



PIERRE VAN MECHELEN a 
C*~) ■ Universiteit Antwerpen 

O 
O 

>>: 




(N 

> 

m 

Results obtained by the HI and ZEUS Collaborations on jet production in photoproduction 
, and deep-inelastic scattering are presented and compared to NLO QCD models. 
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O ■ 1 Introduction 

<j_) ■ The observation of jet production at HERA allows to test the theory of strong interactions, 



quantum chromo dynamics (QCD), and to extract information on the parton content of the 
proton and of real or virtual photons. QCD predicts the production of partons with large trans- 
verse momenta, fragmenting into jets with similar four-momenta. The study of jet observables 
therefore allows to investigate the underlying parton dynamics. Perturbative QCD predictions, 
however, also need the parton content of the proton and photon as input. Jet cross sections can 
thus be used to further constrain the parton density functions (PDF's) obtained by global fits. 

The clustering of final state objects into a few jets is performed by applying a jet finding 
algorithm. The results presented in this paper are all obtained with the so-called inclusive, 
longitudinal invariant k± algorithm^. This algorithm has the advantage of being infrared and 
collinear safe and to be minimally sensitive to fragmentation and underlying event effects. 

Reconstructed data are corrected from detector to hadron level using a full detector simula- 
tion. Next-to-leading order (NLO) QCD predictions are corrected from parton to hadron level. 
This involves the fragmentation of partons into hadrons and secondary interactions between 
partons of the photon and proton remnants. Correction factors are obtained from leading order 
(LO) Monte Carlo (MC) models where NLO effects are modelled by QCD cascades. 

In LO two type of processes are distinguished. In direct processes the exchanged photon 
interacts as a whole with the proton to produce jets. In resolved interactions, the photon is 
treated as a source of partons, one of which produces a hard scattering with the proton, leaving 
a (soft) photon remnant behind. The concept of resolved photons is useful in photoproduction 
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as well as in deep-inelastic scattering (DIS) when Ej, S> Q 2 , where Et is the transverse energy of 
the partons produced in the hard interaction and Q 2 is the virtuality of the exchanged photon. 
To separate direct and resolved enhanced event samples, the momentum fraction of the photon 
entering the hard scattering, be used. On hadron level, the variable xtf , which is 

correlated to the parton level x 7 , is calculated as x 7 et = Yl jets Ej^e^ 3 * 1 /2E^ (with Ej? and 
rf et the jet transverse energy and pseudorapidity and Ej the exchanged photon energy). 

Using QCD factorization, the direct and resolved cross sections for producing N jets, inte- 
grated over phase space, can be expressed as 

dUf j/e (y, Q 2 ) £ f i/p (x p , p? p )a^ N ^\ (1) 

i 

<mf y/e (y, Q 2 ) £ fi/p(x P , f4)fj/v(x>r, ^V^ Njets . (2) 

ij 

Here / 7 / e , fi/ p and /^y 7 are flux factors for photons originating from the electron and for 
partons originating from the proton and the photon, respectively, evaluated at given frac- 
tional momentum-energies and factorization scales. The partonic cross sections 0-7*-*Wjets a nd 
a tj^Njets can ca l cu l a ted in LO or NLO as a function of the strong coupling constant as(fMR), 
where [Ir is the renormalization scale. The choice of renormalization and factorization scales 
will lead to some uncertainty in the predicted cross sections. Different NLO calculations further 
differ mainly in their treatment of infrared and collinear divergences. 

2 Inclusive jet photoproduction 

Figure U shows the first HI data on inclusive jet photoproduction obtained with the k± algo- 
rithm . This analysis comprises two data sets. At large E^, tracker and calorimeter informa- 
tion is emploid to trigger the data acquisition, while at low Ej? a trigger based on the presence 
of the scattered electron in an electron tagger is used, since in this kinematic range the first 
trigger suffers from large backgrounds. The low £jf* data points are corrected for the different 
Q 2 range covered, whereafter they agree well with the high E^f* data in the overlap region. 

The cross section as function of Ej? is displayed in Fig. It falls over more than six orders 
of magnitude from E^ 1 = 5 to 75 GeV. LO QCD predictions underestimate the cross section 
(even with hadronisation corrections), especially at low Ej? 1 , while NLO QCD models reproduce 
the data well after applying hadronisation corrections. Different choices of photon and proton 
PDF's vary at the level of 5 to 10 % and describe the data within errors. 

Figure shows the cross section as a function of r/ je *, in two ranges of E^ 1 . It is ob- 
served that hadronisation corrections increase towards the proton remnant direction. For E^f* > 
12 GeV, the data are well described by NLO prediction, but at lower E^ 1 the data exhibit a 
faster rise towards the proton remnant direction than predicted by NLO QCD. This might be 
explained by a failure of the LO MC models to describe the hadronisation corrections properly, 
by the inadequacy of the photon PDF or by the need for higher order (HO) corrections. 

The ZEUS Collaboration uses the measurement of jets in photoproduction to extract values 
for as (E J T et p. The jet cross section is fitted in each bin of E^ to 

[da/dE' T % = Aas(E^) + B^E^), (3) 

where the constants Ai en B{ are obtained from NLO QCD calculations. A fit of the energy-scale 
dependence of the extracted as(£^ e< ) to the renormalization group equation yields 
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Figure 1: Differential inclusive jet cross sections as a function of (a) and rf et (b). The hatched band shows 
the uncertainty associated with the LAr calorimeter energy scale. The inside and outside shaded bands represent 
the uncertainty from hadronisation corrections and from renormalization and factorization scales, respectively. 



a s (M z ) = 0.1224 ± 0.0001(stat.)l°;°^(syst.)iHo42(^eo.). 
This is in agreement with the current world average of 0.1183 ± 0.0027. 



(4) 



3 Dijet electroproduction 

When the virtuality of the exchanged pho- 
ton is much larger than the natural QCD mass 
scale {Q 2 ^> Aq CD ), resolved photons are in 
principle not needed. However, LO QCD does 
not describe the data on dijet production and 
HO effects therefore need to be included. 

The HI Collaboration compares two ap- 
proaches^. The first approach uses the HER- 
WIG MC model based on DGLAP QCD 
evolution with resolved photons. Both trans- 
versely and longitudinally polarized photons 
are needed to describe the measured cross sec- 
tion, as shown in Fig.[2j The second approach 
uses the CCFM evolution equations as imple- 
mented in the CASCADE MC model© with 
only direct photons. This model is equally 
well able to describe the data. It is therefore 
clear that the ordering of a single DGLAP 
cascade needs to be broken, but whether this 
can be achieved by introducing a second cas- 
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Figure 2: Differential cross section for dijet production 



as a function of a: 7 or y in bins of Q and average ■ 
cade originating from the photon or whether new QCD dynamics is needed can not be decided. 



ZEUS computes the ratio of resolved to direct enhanced components of the dijet cross sec- 
tion^ defined as: 
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c) 150<E 2 < 700 GeV 2 
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Figure 3: The ratio of resolved to direct enhanced components of 
the dijet cross section as function of Q 2 in bins of average E^ 1 . 

4 Inclusive jet electroproduction 
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Figure shows R function of 
Q 2 in bins of average £Sf*. Two 
NLO QCD models are used for com- 
parisons. DISASTER+- 1- 1 ^ only in- 
corporates pointlike photons, while 
JETVIP^ implements both direct 
and resolved photons. Neither de- 
scribes the ratio of resolved to di- 
rect enhanced components well and 
although JETVIP should produce a 
resolved enhanced component it ac- 
tually performs worse than DISAS- 
TER-!- +• ZEUS also observes, how- 
ever, that the direct enhanced com- 
ponent by itself is well reproduced 
by both models. 
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Figure 4: The inclusive jet cross section integrated over 5 < Q 2 < 
100 GeV 2 and 0.2 < y < 0.6. The hatched band around the NLO 
prediction represents the renormalization scale uncertainty. 



Inclusive jet cross sections ob- 
tained by are shown in Fig. 0] 
In the backward region and at all 
r]^l for E J T et > 20 GeV, NLO QCD 
calculations by DISENT^U describe 
the dai et / dEj? 1 distributions well. 
Deviations are, however, visible to- 
wards the proton remnant if both 
E^* and Q 2 are small. These dis- 
crepancies are accompanied by large 
corrections between LO and NLO 
predictions. 

Using a similar method as men- 
tioned before, ZEUS has extracted 
values for as(l!/jf*)^l The result 
is shown in Fig. \5\ and is compared 
to the scale dependence from the re- 
normalization group equation using 
as input value the result obtained 
from a fit to the measured da^ et /dQ 2 
distribution for Q 2 > 500 GeV 2 : 



a s (M z ) = 0.1212 ± 0.0017(stat.; 
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5 Conclusion 
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Figure 5: The asiE^ 1 ) values as determined from a 
QCD fit of the measured jet cross section. 



The HI and ZEUS Collaborations have 
measured jet production in ep collisions with 
real and virtual photons in a large kinematic 
range. Cross sections are obtained with high 
accuracy and fall over more than six orders of 
magnitude ctS Si function of transverse energy. 

Competitive values for the strong coupling 
constant as are obtained which are in agree- 
ment with the current world average. 

NLO QCD calculations do a very good 
job in describing jet cross section, with excep- 
tions for forward jets at low Q 2 and Ej? and 
for the ratio of direct to resolved enhanced 
components in dijet production. 

This wealth of new high-precision data 
can, and should, be used in global fits of pho- 
ton and PDF's. 
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